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Abstract - The relationship between cirrus particle size and 
cloud top temperature is surveyed on a near-global scale. 
The cirrus particle size is retrieved assuming ice crystals are 
hexagonal columns and the cloud top temperature and the 
radiaces in channel 1 and 3 of AVHRR used to retrieve ice 
particle sizes are from ISCCP product. The results show that 
for thick clouds over North America, the relation between 
cirrus particle size and cloud top temperature is consistent 
with a summary of this relationship based on aircraft 
measurement over that region for thick clouds. However, 
this relationship is not universal for other regions especially 
for tropical zone, which has been found by other in situ 
measurements. 
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parameterize the ice crystal size as Heymsfield and Platt [6] 
suggested. While these data supply valuable information 
about this relationship, they are all regional measurements 
during a short time period for mostly thick cirrus clouds. 


This study investigates this relationship based on the 
cirrus particle size data retrieved from ISCCP CX data. We 
found that for thick clouds (t > 10), similar to those found by 
aircraft measurements, most of the regions over the globe 
show positive relationships between cloud temperature and 
cirrus ice crystal sizes. However, if relatively thin clouds are 
included ((t >3), this correlation becomes negative for 
tropical areas but remains mostly positive for midlatitudes. 

RETRIEVAL OF ICE CRYSTAL PARTICLE SIZE 


INTRODUCTION 

Cloud microphysical parameterizations have attracted a 
great deal of attention in recent years due to their effect on 
cloud radiative properties [1] and cloud-related hydrological 
processes in large-scale models [2] The parameterization of 
cirrus particle size has been demonstrated as an indispens- 
able component in the climate feedback analysis [3]. In 
parameterization schemes, relationships between cloud 
temperatures/optical thickness and cirrus ice crystal sizes are 
of critical importance because the microphysics of cirrus 
clouds modifies the relationship between cloud optical depth 
and cloud ice/liquid water path. Climate models not 
accounting for this relationship cannot correctly predict the 
temperature dependence of infrared emittance of cirrus 
clouds [4] and thus arrive at wrong conclusions about 
climate change. Many efforts have been made to investigate 
the cloud temperature-cirrus crystal size relationship and 
controversial conclusions have been reached. For example, 
twenty flights of aircraft measurements over North America 
during winter times [5,6] show that the cloud-particle mode 
radius increases with temperature. Similar results were 
found by in situ measurement during CEPEX at tropics [7] 
and by surface lidar observations [8]. However, 
investigations using radar measurements [9], Doppler radar 

and IR radiometry [10], and millimeter- wave radar [11] 
found a negative correlation between ice crystal size and 
temperature and conclude that temperature cannot be used to 


The detailed methodology of retrieving ice crystal size is 
presented in another paper [12]. In that paper, the method 
for retrieving cirrus particle size information on a near- 
global scale (50°S to 50°N) using currently available satellite 
data from ISCCP is described. To retrieve cirrus particle 
size, we use a radiative transfer model that includes all major 
absorbing gases and cloud scattering / absorption to compute 
synthetic radiances as a function of satellite viewing 
geometry. Ice clouds are determined by cloud top 
temperatue T c £ 240K. Shapes and orientations of ice 

crystals are assumed to be hexagonal columns and plates 
randomly oriented in the atmosphere. Although ice crystals 
may have complicated shapes such as bullets, rosettes, 
aggregates, irregular and quasi-spherical particles and show 
preferred orientations, it is difficult to determine which 
shape and/or orientation is dominant in a specific clouds by 
current remote sensing instruments. On one hand, we 
realize the possible effect of non-hexagonal particles and 
preferred orientations. On the other hand, we also notice 
that the difference of phase functions between hexagon and 
other irregular shapes is much smaller than that of non- 
spherical and spherical particles [13]. Therefore, we adopt 
the phase functions for hexagonal columns in the retrieval 
and estimate the possible effect of non-hexagonal and 
preferred orientation by comparing the retrieved particle 
sizes with in situ measurements. Ray tracing techniques are 
used to calculate phase functions for different size 
distributions. Five different size distributions from 


Acknowledgement . We wish to thank Dr. Ou and Takano for supplying phase functions used in this study. This research was supported by NASA Grant Nos. 
NAGW-3922, NAGW-3788; was partially funded by DOE through the NIGEC Great Plains Regional Center at the University of Nebraska- Lincoln (DOE 
Cooperative Agreement No. DEFC03-90ER61010). Financial support does not constitute an endorsement by DOE of the views expressed m this paper. This 
research was also supported by the NASA Climate Program managed by Dr. Bob Curran. 


0-7803-3836-7/97/S1 0.00 © 1997 IEEE 


1760 











cirrus (D e =47.6 jim), -40°C cirrus (D e =64.1 jam), Nov. 1 
cirrus (D e =75.1 jam) and Cirrus uncinus (D c =123.6 pm). 
Phase functions of these five size distributions for channels 1 
and 3 of AVHRR are applied in the model for calculations of 
multiple scattering.The model results have been validated 
against clear sky observations and are consistent with the 
observed radiance range under cloudy conditions. 

The results of this near-global survey show that the peak 
frequency of cirrus crystal effective diameter D e for hexagons 
is about 60 pm, the mean value corresponding to a 30 pm 
radius for an equivalent sphere. The survey also reveals that 
about 15% of the small crystal sizes are less than 23 pm in 
effective diameter. The geographical distribution of small 
crystal sizes is consistent with the results of Prabhakara et al. 
[14]. There are no significant differences between particle 
sizes of continental and maritime cirrus. 

RELATION BETWEEN ICE CRYSTAL SIZE 
AND TEMPERATURE 

Figure 1 shows the distribution of correlation 
coefficients between ice cloud crystal size and cloud top 
temperature for thick clouds (t > 10). Ice clouds over most 
of the regions show positive relation between ice crystal size 
and cloud top temperature, consistent with the results of 
Heymsfield and Platt (1984). However, this relation changes 
when thin clouds are included into the statistics. Figure 2 
shows that for ice clouds with (x £ 3, most ice clouds over 
tropical and sub-tropical regions show negative relation 
between ice particle sizes and cloud top temperature. The 
possible explanation includes different cirrus systems over 
tropics and midlatitudes. The strong convective motion in 
tropics brought large ice crystal particles up to very high and 
cold regions and spread out to form vast and relatively thin 
anvils. Further information about the underlying mechanism 
and possible explanations are under investigation. 
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